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Abstract: We evaluate all phase space master integrals which are required for the total 
cross section of generic 2—7-1 processes at NNLO as a series expansion in the dimensional 
regulator e. Away from the limit of threshold production, our expansion includes one order 
higher than what has been available in the literature. At threshold, we provide expres- 
sions which are valid to all orders in terms of F functions and hypergeometric functions. 
These results are a necessary ingredient for the renormalization and mass factorization of 
singularities in 2 — )■ 1 inclusive cross sections at N^LO in QCD. 
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1. Introduction 



Signals of novel physics, such as the recently discovered Higgs boson by the ATLAS [1] and 
CMS [2] collaborations, as well as standard candles, such as Drell-Yan production of elec- 

troweak gauge bosons [3, 4], have been at the epicenter of experiments at hadron colliders. 
The analysis and interpretation of data relies crucially on high precision theoretical esti- 
mates, based on perturbation theory, for the rates of 2 ^ 1 processes. Most importantly, 
higher order corrections due to Quantum Chromodynamics (QCD) are generally large at 
LHC energies, altering the central value and lowering the theoretical uncertainty of many 
observables. 

The state of the art in fixed-order perturbative QCD for 2—7-1 processes is next-to- 
next-to-lcading-order (NNLO). The first inclusive cross section computation at this order 
was performed more than two decades ago for Drell-Yan production [5]. This calculation 
was followed by several computations for inclusive Higgs production cross sections [6, 7, 8, 
9, 10, 11, 12, 13, 14, 15] with the emergence of powerful computational techniques [6, 7]. 

The precision of NNLO cross sections is quite high as indicated by the study scale 
variations. For example, the uncertainty due to scale variations in Higgs production via 
gluon fusion is less than 10% and in Drell-Yan production less than 2% (sec, for example, 
ref. [16] and ref. [17]). However, a more reliable estimate of the uncertainty due to the 
truncation of perturbation theory would be the evaluation of yet another term in the 
series expansion. Now that a Higgs boson has been discovered, measuring its properties 
and couplings precisely becomes a priority. With this prospect of precision physics in 
the Higgs sector and the paramount importance of Higgs coupling measurements for our 
understanding of physics at high energies, it is evident that many important observables 
will need to be validated at an even higher order in perturbation theory, i.e. at next-to- 
next-to-next-to-leading order (NNNLO or N^LO). 

In this paper, we present one of the ingredients needed for 2 — )• 1 inclusive cross sections 
at N^LO. In particular, we compute the double-real and real-virtual master integrals for a 
generic 2—7-1 scattering amplitude at one order higher in the dimensional regularization 
parameter e = (4 — D)/2, with D the space-time dimension, than previously known [7]. 
The paper is organized as follows. In section 2, we motivate our computation by explaining 
where they enter the calculation of a N^LO observable. In section 3 we provide some details 
about the method we use. Sections 4 and 5 contain the results and we conclude in section 
6. 

2. NNLO contributions to N^LO observables 

In this paper we are concerned with the computation of the fully inclusive cross section 
for the production of a massive state X. In perturbation theory this cross section can 
schematically be written as 



oo 
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(2.1) 
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where 

T = and z = (2.2) 

o s 

with S being the total center of mass energy, s the partonic center of mass energy carried by 
the momenta of the incoming partons, pi and p2, and mx the (on-shell) mass of particle X. 
The bare strong coupling constant is denoted as a^s'^ , ff^ are the bare parton distribution 
functions (PDF's) and aij^n is the partonic cross section for X + n jets, which in turn 
admits a perturbative expansion in the number of loops. may be written as a phase 
space integral over the squared X + n jets amplitude as follows, 

where we define the measure of the phase space volume for a massive particle X of mass 
mx and n — 1 massless particles by 



where 




(2.5) 



Using this notation, the Lorentz invariants appearing in the phase space integrals are 
defined through 

Sh.-in = iPil..inf ■ (2-6) 

While the inclusive cross section ax must obviously be finite order by order in pertur- 
bation theory, the individual pieces contributing to a given loop order are divergent. Final 
state infrared (IR) divergences cancel mutually between the real and virtual corrections, 
whereas ultraviolet (UV) and initial-state IR divergences have to be dealt with by replacing 
the bare coupling and PDF's by their renormalized counterparts, which requires the in- 
troduction of explicit counterterms proportional to poles in the dimensional regularization 
parameter e multiplying lower order coefficients of the cross sections. As an example, if 0"^^^ 
denotes the i-th order correction to the cross section, the N^LO UV/PDF counterterms 
can schematically be written as 

<5a(3) ~ X + (^Ca + ^Cs) x a^'^ + (^^C, + ^C, + icg) x , (2.7) 

where the x may indicate a convolution in the case of the PDF counterterms. As a 
consequence, the poles of the counterterms produce finite contributions to the N^LO cross 
section from the higher orders in the e expansion of a^^\ 

The aim of this paper is to provide the NNLO master integrals which are required 
in (T*^^) to one order higher in the e expansion. At this order in perturbation theory three 
different contributions need to be taken into account. As an example, the NNLO correction 
to the partonic cross section for gg ^ H gets contributions from the following three types 
of interference diagrams, 
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• double-virtual: the two-loop production amplitude interfered with the Born am- 
plitude as well as the square of the one- loop amplitude for gg H, e.g., 




• real-virtual: the one-loop production amplitudes for gg — t- Hg, gq — ?• Hq and 

gq —7- Hq interfered with the corresponding Born amplitudes, e.g.. 




• double-real: the square of the Born amplitudes for gg — t- Hgg, gg — )■ Hqq, gq — >• 
Hgq, gq Hgq, qq Hqq, and qq Hqq, e.g.. 




First, the double-virtual corrections decouple from the phase space integration, reducing 
the computation of the two-loop phase space integrals effectively to the computation of 
the two-loop QCD form factor. The master integrals for the two-loop QCD form factor 
have been evaluated in ref. [18, 19, 20] to all orders in e in terms of T functions and 
hypergeometric functions, which can easily be expanded in e using the HypExp package [21]. 
For this reason we will not consider these integrals any further, but we simply list them 
in appendix A.l for completeness. Next, in ref. [7] it was shown how to compute the 
phase space integrals describing the real-virtual and double-real corrections as a Laurent 
expansion in the dimensional regulator up to order e'^. More precisely, it was shown that 
all the relevant integrals can be reduced to a small set of master integrals, which were 
computed to the order in e required for NNLO computations. In the rest of this paper we 
compute all the master integrals to one order higher in the e expansion, thus preparing the 
ground for using them in computations beyond NNLO. 

3. Phase space integrals from differential equations 

In this section we give a brief account on how to compute the master integrals for the real- 
virtual and double-real phase space integrals. Our computations are based on the method 
introduced in ref. [7], which we shortly review in the following. In ref. [7] all real- virtual 
and double-real topologies were expressed as a two-loop forward scattering amplitude with, 
respectively, two and three cut propagators. Such cut integrals are related to the desired 
phase space integrals via Cutkosky's rule. 

Disc ^5 ^ = 27Ti6+(p'^ - ra^) = Imdip^ - n?) d(f) . (3.1) 

pA _ jj^Z _|_ 
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Eq. (3.1) allows us to interpret a phase space integral as a loop integral, and thus we can 
apply results developed for the computation of Feynman integrals. In particular, these 
two-loop integrals can be reduced to a set of master integrals by virtue of the Laporta 
algorithm [22], using integration-by-parts (IBP) and Lorentz invar iance identities, as im- 
plemented in the program AIR [23] , with the additional constraint that integrals with some 
of the cut propagators not present vanish and need not to be considered in the reduction 
process [7]. 

The master integrals themselves are calculated using the method of differential equa- 
tions [24, 25]: we differentiate each master integral under the integration sign with respect 
to a kinematic invariant, in our case the mass or equivalently the variable z. 

During the differentiation process, loop integrals with higher powers of propagators 
are produced. They can again be expressed in terms of the master integrals themselves via 
the IBP identities and the system of differential equations closes upon itself. If we denote 
the master integrals by Fj(2;, e), we arrive at a set of first order ordinary linear differential 
equations for the master integrals, 

— F,(^, 6) = ^ (z, e) F,(z, e) , (3.2) 

i 

where the Cjj(z, e) are rational functions of z and e which have poles in z at most at z = 
and/or z = ±1. In the case where the set of differential equations is triangular order by 
order in e, the set of equations may be solved order by order using standard techniques 
for ordinary differential equations. In particular, if oj{z) is a homogeneous solution to the 
linear differential equation 

-y{z) = A{z)y{z) + Biz), (3.3) 
then the solution to the inhomogeneous differential equation is given by the integral 

yiz) = a uiz) +u{z)j dz' ^ , (3.4) 

where a G M is an integration constant. The poles appearing inside the integral in eq. (3.4) 
are determined by the poles in the coefficient A{z). Thus in our case the solution for the 
master integrals can be expressed through (iterated) integrals with poles at most at ^; = 
and z = ±1. Integrals of this form lead naturally to harmonic polylogarithms [26], defined 
recursively through the iterated integrals 

H{au. 

where G {—1, 0, 1} and 

f{-ht) = 

In the case where all the are zero, the integral (3.5) is divergent, and we define instead 

H{On;z) = -Aog^z. (3.7) 
n! 



i;z)=[ dtf{ai,t)H{a2,...,an;t), 
Jo 



(3.5) 



(3.6) 
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In this way we can express all the master integrals, up to the desired order in e, in terms of 
harmonic polylogarithms up to weight four, which can be evaluated numerically in a fast 
and accurate way [27, 28, 29, 30, 31]. 

In order to fully determine the solutions of the differential equations, we need to fix 
the integration constants. This can be achieved by requiring the master integrals to take 
particular values, computed separately by other means, for some special value of z. Since 
< z < 1, there are two natural choices for such special values. The case z = physically 
corresponds to the situation where the produced particle X is massless. However, the limit 
z — > is in general not smooth, as the amplitude might develop new infrared poles in this 
limit. For this reason we choose the point z = 1 as an initial condition to the differential 
equation, which corresponds to the soft limit in which the momenta of all the final-state 
partons vanish. As this limit is important not only as an initial condition for the differential 
equations, we review the structure of the phase space integrals in the soft limit in the next 
section. 



3.1 Soft limits of master integrals 

In the previous section we argued that the point z = 1 is a good initial condition for the 

differential equations satisfied by the master integrals. There is however another reason 
why the point z = \ deserves special attention. Indeed, some of the master integrals contain 
simple poles at z = 1, and without loss of generality we can always write such integrals in 
the form 



nezx 



(l-z) 



1+ne) 



where n is an integer and the function Fn{z,e) is finite at z 
divergence associated to the soft singularity by expanding 



(1 



-l+ne 



ne 



k=o 



where 



Vk{x) 



log^(x) 



and the + indicates the common plus-prescription, 

log^Xx) 



dxVkix) fix) = C dx \f{x) - f{0) 

Jo Jo ^ 



(3.8) 

We can isolate the 
(3.9) 

(3.10) 
(3.11) 



It then becomes apparent that, if we want to know F{z, e) up to a certain order in the e 
expansion, we need to know -F„(l, e) at one order higher in the expansion in e than Fn{z, e), 
since Fn{l,e) is multiplied by the pole in eq. (3.9). More generally, even if a given master 
integral has no pole as z approaches 1, when appearing inside an amplitude it might get 
multiplied by a coefficient such that the product develops at most a simple pole at z = 1. 
Hence, we have to use eq. (3.9) to perform the e expansion and therefore, by the same logic, 
we always need to know the values of the master integrals at z = 1 to one order higher in 
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the e expansion. Thus, there is an independent physics motivation to study the soft hmits 
of the master integrals in some detail. 

We define the soft limit of a master integral F(z, e) as the (unique) function 

Fnie) 



F^(z,e) = ^ aeZ>i, (3.12) 

such that 



Mm ' = 1 . (3.13) 



We stress that the right-hand side of eq. (3.13) is valid to all orders in e. In addition we 
define the hard part F^{z,e) of a master integral by 

F(z,e) = F^(z,e)+F^(z,e). (3.14) 

In sections 4 and 5 we will determine the soft limits of the master integrals directly to all 
orders in e, without the help of the differential equations. Our knowledge of the soft limits 
combined with eq. (3.13) as an initial condition for the differential equations then allow us, 
at least in principle, to compute the analytic expressions for the master integrals up to any 
order in e. In the rest of this paper we perform this task explicitly for all master integrals 
up to transcendental weight 4. 

4. Double-real master integrals 
4.1 Definitions and conventions 

In this section we compute the analytic expressions for the master integrals for the double- 
real emission contributions. In ref. [7] it was shown that all phase space integrals for the 
production of a massive particle and two massless particles can be reduced to a linear 
combination of 18 master integrals, which we choose as follows: 



Pi 



P2 




yPi r 

^ = / d$3 = s'-'^7'(e)Xi(z,e), 

^P2 J 



(4.1) 



P2^ 



*^P2 



Pi 



s-i~2^p(e)X2(z,e) 



■S234S134 



(4.2) 



Pi 



P2 




Pi 



P2 



S14S23 



P(e)X3(z,e), (4.3) 
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P2^ 



*-Pl 



d^3 



S234'523Sl3'Sl34 



(4.4) 



P2^ 




*-P2 



d^3 



s--'-^^r{e)X,{z,e), 



S234S24Sl3'Sl34 



Pi 



(4.5) 



P2^ 



*-P2 



*-Pl 



S13S23S14S24 



s-'-'^r{e)Xe{z,e), 



(4.6) 



Pi 

P2 




Pi r 

P2 J ■S123S124 



'^^3 -1-2. 



s-'-'^r{e)Xr{z,e), (4.7) 



P2^ 



*-P2 



*-Pl 



d^3 



s-^-^^V{e)Xs{z,e), (4.8) 



Sl23'Sl24'Sl4Sl3 



Pl^ \ 71*- Pi 




d^3 „-l-2e 



S14S23 



s''~'^r{e)X,{z,e), 



P2^ \^P2 



(4.9) 



Pl^ 



P2^ 



\ 




\ 






\ 




\ 



*-pl 



s-3-2^P(e)Xio(z,e), (4.10) 



Sl23'Sl24'Sl4'S23 



*-P2 



Pi 
P2 




Pi 
P2 



J d$3S34 = s'-'^ Vie) Xii(z, e) , (4.11) 
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P2^ 



^P2 



■S234'Sl23 



(4.12) 




(4.13) 



(4.14) 



(4.15) 



— \ V t> Pi 



P2^ ^\^> P2 



I ^ ^_3„2, ^^^^^^ ^ 



Pl^ 



P2> 



\ 




\ 






\ 




\ 



*-P2 



*-Pl 



-3-2e 



S34'S13'S234S123 



P(e)Xi7(z,e) 



(4.17) 




P(e)Xi8(z,e), (4.18) 



The solid lines denote the massive particle and dotted lines represent numerator factors. 
In all integrals we have pulled out an overall dimensionful scaling factor as well as an overall 
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normalization 

Note that the integral Xi is just the phase space volume, which can be written to all orders 
in e as 

Xi(^, e) = (1 - zf-^^I^^^,F,{l - 6, 2 - 2e; 4 - 4e; 1 - ^) . (4.20) 

For the other integrals no similar all-orders formulae are known to us. We therefore compute 
these integrals using the differential equation technique. It turns out that the system of 
differential equations for the double-real master integrals are triangular order by order in 
e, except for Xi and Xn, which can be decoupled order by order by replacing Xn by the 
linear combination 

Xii„ = ^[(l + 2)Xi-Xn] . (4.21) 

After this change of variables the system is triangular order by order in e, and we can 
immediately solve for the master integrals, provided that we know the soft limits of the 
master integrals which serve as an initial condition. The computation of the soft limits will 
be the topic of the rest of this section. 

4.2 Soft limits of double-real master integrals 

To evaluate the soft limits of the double-real master integrals we use the following exact 
parametrization for the 2^3 phase space, 

s'^z^K^ sin^(7rx2) JJ yi 1 , 



4 

SZ K 



X 



K, 



i=l / 

with z = 1 — z. In this parametrization the propagators of the massless partons read 

Si3 = -SZKyi, S23 = -SZKy2, 

su = -szKys, S24 = -SZKyi, (4.23) 

where 

i = z (yiy4 + ^2^3 + 2 cos{x2■K)^Jylymy) , (4-24) 

and 

'^ = ^— e[l,2). (4.25) 

This parametrization may be derived from the "energies and angles parametrization" of 
ref. [32] via the following transformation 

yi = X1X3, y2 = X1X3, y3 = xiXi, ?/4 = X1X4, , (4.26) 
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with Xi = 1 — Xi- In addition, the master integrals depend on the following denominators 

Sl34 = Sl3 + Si4 + S34 , S234 = S23 + -524 + ^34 , 

S123 = S12 + S23 + Sl3 , Sl24 = S12 + S24 + Sl4 . 

We now take the soft limit, z ^ 1. It is easy to see that in the soft limit we have the 
relations 

lim si23 = lim S124 = 512 , lim k = 1 , 

lim S134 = si3 + si4 , Hm S234 = S23 + ^24 ■ 

Z-^l 

If we now perform the change of variables 

yi = hts, y2 = tit4, ys = hh, 2/4 = hh ■ (4.29) 

It is then easy to see from eqs. (4.23) and (4.29) that in the soft limit all the invariants, 
except for S34, take a fully factorized form. In particular the terms S13 + S14 and S23 + ^24 
effectively reduce to ti and ti respectively. This construction therefore allows one to derive 
all the soft limits of almost all master integrals, with the exception of X17 and Xig, in 
terms of simple Beta-functions, 

B(x,.) = /'*t-(l-«r' = ™. (4.30) 

The integral X17 may be trivialized in a similar way in the "hierarchical parametrization" 
of ref. [32]. Unfortunately, we are not aware of any parametrization which allows one to 
factorize all denominators of Xig simultaneously. We proceed by expressing the integral in 
terms of energies and angles, i.e., by writing pi = £'j(l, nj), where Hi is a unit vector in the 
direction of the spatial component of pi. The integral over the energies of ps and p4 yields 
just another Beta-function, and we arrive at 

P(6)X,3 - 4(47r) (i_,)i+4. J (l + cosM(l-cosM(l-cosM ' ^ ' 

where dil^'^ is differential volume element of the d — 1 dimensional solid angle of pi and 
the angles are defined as cos dij = UiMj. We first use the result [33] 



(d-i) 
4 



cos 014)(1 — cos ^34) 

^ j^(d-3) r dgi3(singi3)'^-3rf034(sin034)'^-" .4 32. 

7o (1 — COS0i4)(l — COS014COS013 — COS(?!)34Sin014Sin^l3) 

~ 7(47r) — ——2^1 1,1,1 -e, 



4' ' er(-2e) 
The integral then becomes 

n^)^i, - -(i_,)i+4. ,r(-4e) ^ 



y X (4-33) 

-dMsin%y^ P 1 1 l + cos^is' 

— — T 2i'i 1, 1; 1 - e; — 

1 + cos 6'i3 V 2 
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Changing variables to y = l+S^Aa^ gq (^4.33^ may be written as 

neMs = (^^^P. '^^''^'i^ / dyy-'-'i^ - y^.F, (1, 1; 1 - y) 
1 (47r)-3+2^ r(-e)2 



2 (1 - e2r(-4e) 



F2(l,l,-e;l-e,l-2e;l), (4.34) 



where wc used the recursive definition of the pFq function. Hence we arrive at expressions 
valid to all orders in e for the soft limits of all the double-real master integrals. We observe 
that in all cases, except for Xfg, the results are proportional to the soft limit of the phase 
space volume, 

Xf(z,e) = (l-.)3-4^^|^, (4.35) 
the constant of proportionality being a rational function of z and e. We therefore define 

Xf{z,e) = Si{z,e)Xf{z,e). (4.36) 
In this normalization the results for the soft limits of the master integrals read 

Si{z, e) = Sjiz, e) = Sna{z, e) = l, (4.37) 



2(3 - 4e) 
(l-26)(l-z) 

,(l -26)(3- 46) 

£2(1 - z) 

(1 - 2 e) (3 - 4e) (1 - 4e) 
£3(1 - zY 



S2iz,e)= "1^^^;^^ (4.38) 



Ssiz, e) = Ssiz, e) = Sg{z, e) = Sio(z, e) = 2 ^ ^ , (4.39) 



S4(., e) = S,{z, e) = -2 ^^"^^^^r 2 " ^ > (4-40) 



Se(.,e)^-8 (^-^-)g-^y-^-) , (4.41) 
e'*(l — z) 

Si2(^,e) = ^, (4.42) 

"^3(-'-) = - (l-l)('l-.) ' (^-^^^ 

Si4(^,6) = ^^^, (4.44) 



(l - 26)(3 - 46) 
e2(l - z) 

(1 - 2 6) (3 - 46) (1 - 46) 

63(1 - Z) 



s,,{z,e)= '^ (4.45) 



Si6(., e) = 'X ' ' (4-46) 



S.(.e)^-2ii^M|^i^, (4.47) 
e'^(l — z) 
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18 132 1 

- ^ + — + - (I2C2 - 288) - 88C2 + 6OC3 + 192 



(1 - zY 



'S72 

+ ( — Cl + 192C2 - 44OC3 ) e 



+ - 



+ O(e^) 



2728 



-Cl + 24C3C2 - I28C2 + 444C5 + 960C3 e 



(4.48) 



One might wonder whether our result for Sis, eq. (4.48), can be expressed to all orders 
in e through ratios of T functions only. Without loss of generality we can assume that in 
that case the 3F2 function entering Sig could be written in the form, 



Rie) l[T{l±nk€) 



mk,nk e 



(4.49) 



where R{e) is a rational function. In the following we argue that the form (4.49) can be 
excluded for Sig, i.e., that 3F2 (1, 1, — e; 1 — e, 1 — 2e; 1) cannot be reduced to a ratio of T 
functions only. Indeed, if we assume that the 3F2 function in eq. (4.48) can be written in 
the form (4.49), then, because of the well-known relation 



r(l + ne) = exp 



k=2 



-ne) 



k 



(4.50) 



its e expansion will to all orders in e only involve zcta values of depth one, i.e., it will 
be possible to write its e expansion without multiple zeta values (MZVs). Multiple zeta 
values are known to satisfy many identities among themselves, and it is conjectured that 
all the identities among MZVs are generated by the double shuffle relations. Using the 
double-shufHe relations, it can be shown that the first time a MZV cannot be written as a 
polynomial of MZVs of depth one is at weight 8, and this irreducible MZV can be chosen to 
be ^5,3- Using HypExp we can perform the expansion of 3F2 (1, 1, — e; 1 — e, 1 — 2e; 1) up to 
0(e^). The result obtained by HypExp is a complicated combination of MZVs of weight 8. 
Using the PSLQ algorithm we fit this combination of MZVs to a basis of MZVs of weight 
8, which we choose to be C2 Q, Cs Csj Cs.a}- We obtain 



3^2(1, l,-e;l-e,l-2e;l)|^8 



7489 
'I25 



74C3C5-yC5,3. (4.51) 



We see that at 0{e^) the expansion of 3F2 (1, 1, — e; 1 — e, 1 — 2e; 1) involves an irreducible 
MZV of weight 8. Conversely, we checked that if we remove this MZV from our basis at 
weight 8, the PSLQ algorithm fails to converge. We thus conclude that the 3F2 function 
appearing in Sig cannot be written as a ratio of F functions of the form (4.49). 

Having at our disposal both the differential equations and the soft limits for all the 
master integrals, we can in principle solve the differential equations to any order in e. 
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We have explicitly performed this task and computed the e expansion of all the master 
integrals up to transcendental weight 4. As the results are rather lengthy, we do not show 
them here explicitly but collect them in appendix A. 3. The results only involve harmonic 
poly logarithms up to weight 4. In ref. [34] it was shown that up to weight 4 almost all 
harmonic polylogarithms can be expressed through classical poly logarithms only, and that 
only three new functions are needed. In our case only one new function appears, which we 
choose to correspond to the harmonic polylogarithm 

Li2,2(-1, z) = -H{0, 1, 0, -1; z) . (4.52) 

We have checked that our results agree with the results of ref. [7] up to terms in the e 
expansion with transcendental weight at most three. In addition we have checked all our 
results numerically for a few nontrivial values of z. For this purpose we used the method 
of non-linear mappings introduced in ref. [32] . 



5. Real- virtual master integrals 
5.1 Master integral definitions 

In ref. [7] six master integrals for the real-virtual contributions to a 2 —t- 1 partonic cross 
section were identified. Defining the prefactor 

l(47r)- r(l + 6)r(l-e) 
= 2^ r(2 - 2e) ' ^'-'^ 



as well as the one-loop integrals 

Bub(p2) : 



d^k 1 

i7r^/2 A;2(A; + p)2' 



d^k 1 

Tri(s23,Sl23) = / — nrtTiyr-^ ^fZTl A2' (^-2) 



Box(si2,S23,Sl23) 



k^^k+piY{k+pi2?{k+pi2zf 



ITT 



the real-virtual master integrals may be defined as the following two-body phase space 
integrals. 



P2^ 



-Pi 



-P2 



I 



c?$2Bub(s 



13; 



s-2^Q(e)Yi(z,e), 



(5.3) 




J d^2Buh{m\) = s-2^Q(e)Y2(z,e), (5.4) 
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P2' 



I 



(5.5) 




j dci)2Tri(si3,m2,) = s-i-2^Q(e)Y3(z,e), 



(5.6) 



P2^ 



*-P2 



S23 



*-Pl 



(5.7) 




j d$2Box(si3,S23,"ix) = s"^"^'Q(e)Y6(z,e) 



(5.8) 



5.2 Soft limits of real-virtual master integrals 

For the evaluation of the soft limits of the real-virtual master integrals we employ the 
strategy presented in ref. [35]. This strategy involves the following steps: 

1. Substitute known expressions for the one- loop master integrals which are valid to all 
orders in e. 

2. Apply analytic continuation formulae to the hypergeometric functions which express 
the box integrals, if required. 

3. Evaluate the integrals at z = 1. 

We find that the soft limits of all master integrals can then be trivially evaluated as Beta- 
functions: 



' ^ ^ ^ er(2-3e) 



Y^{z,e) = (1 - z)^~2e cos(7re) 



r(i-6)^ 

er(2 - 2e) 



Yf(z,e)=Yf(z,e), 



(5.9) 

(5.10) 
(5.11) 
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(5.12) 



t(z,e) = -3(l-z) 



i_3. r(2-2e)r(l-e) 
e3r(l-3e) 



(5.13) 




{z,e) = -2{l-z) 



-l-4e 



cos(7re) 



r(l - 2e)r(2 - 2e)r(l - e)r(l + e) 
e3r(l -4e) 



(5.14) 



Note that only the real part of the master integrals is shown, as the imaginary part never 
enters the computation of a physical observable. 

As in the double-real case, we list the results for generic z that we obtained by solving 
the differential equations in appendix A. 2. The expressions agree with the results of ref. [7] 
up to terms in the e expansion with transcendental weight at most three. The new orders 
were checked numerically using the approach of [35]. This strategy is as follows: 

1. Substitute known expressions for the one-loop master integrals which are valid to 
all-orders in e. 

2. Apply analytic continuation formulae to the hypergeometric functions which express 
the box integrals, if required. 

3. Expand the hypergeometric in terms of polylogarithmic functions. 

4. Expand the real emission singularities (these are factorized here) in terms of delta- 
and plus-distributions. 

5. Evaluate the coefficients numerically. 
6. Conclusions 

We have given results for all double-real and real-virtual master integrals needed in the 
calculation of a generic 2—7-1 inclusive scattering cross section. The analytic results for 
the master integrals are provided in Maple input form as ancillary material to the arXiv 
submission. The previously known expressions have been supplemented with one more 
order in the e expansion. Together with the previously published double-virtual results [20] 
listed in appendix A.l, all contributions to a N^LO cross section involving NNLO master 
integrals can thus be computed, as outlined in section 2. 

We have divided the master integrals into a soft and a hard part, each of which is given 
separately. This separation is motivated in a twofold way. First, the soft limits provide 
a natural boundary condition to solve the differential equations for the master integrals. 
Second, divergences due to soft radiation have to be subtracted, leaving poles at z = 1 and 
thus requiring the knowledge of the master integral to one order higher in the e expansion, 
as explained in section 3.1. 

We have found the method of differential equations to be very adequate for this type 
of task, and see no obvious problems in applying the same techniques to integrals of higher 
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complexity. Such integrals could include the triple-real master integrals which constitute 
another part of a 2 — )• 1 N^LO cross section. In that case, the main difficulty might be 
the vastly increased amount of master integrals, as well as a much more intertwined set 
of differential equations due to more complicated IBP identities for the triple-real topolo- 
gies. Once these difficulties overcome, one should be able to perform analytically all the 
phase space integrals appearing in 2 — )• 1 cross sections at N^LO, and thus, by combining 
them with the recently computed three-loop QCD form factor [36, 37], arrive at the first 
prediction for an LHC observable computed perturbatively up to N^LO. 
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A. Results for the master integrals 

In this appendix we collect the analytic results for the hard part of the real-virtual and 
double-real master integrals omitted throughout the main text. As the imaginary part of 
the master integrals does not enter the computation of any physical observable, we only 
show the real part. For completeness we also include the analytic expressions for the two- 
loop virtual master integrals, which can easily be obtained from the expressions for the 
two-loop form factor master integrals of ref. [18, 19, 20]. 

A.l Results for the double- virtual master integrals 




(A.l) 




(A.2) 




(A.3) 
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P2^ 



5(1 — z) cos(27re) 



r(l - 2e)r(l - e)^r(l + e)r(l + 2e) 
2(1 - 2e)e2r(2 - 3e) 



(A.4) 



">- 



= 5{l-z) cos(2^e) r(l + ef[\-l-C2--C^- ^Cl (A.5) 



P2 



g4 g2' 



+ 



IO2C2C3 - 117C5 



A. 2 Results for the real-virtual master integrals 



Yf(z,6) 
Yf(z,6) 



Yf(z,e) 



0, 

(z--l)Yf(z,6) 
0, 



>l-2e 



cos(7re) 



r(i 



er(2 - 2e) 



(A.6) 
(A.7) 
(A.8) 



2e . 



Li2(2) H C2 



+ -5Li3(l -z)- 4Li3(z) 



+ 4Li2(z) log z - 5Li2(z) log(l - z) + 2C2 log z + 6(2 log(l - z) + 



log^ z 



log2(l - z) log z + 4C3 - 2Li2(z) - log2 z + 2C2 + e 2C2Li2(z) 

19, 



+ 14Li4 ( 1 - - 



(A.9) 



5Li4(l -z)- AUi{z) + 5Li2(z) log2 z + y Li2(z) log2(l - z) 



- 14Li2(2) log(l - z) log z - 14Li3(l - z) log z - 10Li3(z) log z 

19 

+ 19Li3(l - z) log(l -z) + 14Li3(z) log(l -z) + 7C2 log^ z - —Q2 log'(l - z) 



+ 10C3logz-14C3log(l-z) + 



71ogS 7 



log(l — z) log^ z 



19 7 2 

+ — log3(l -z)\ogz-- log2(l - z) log2 z - -Cl + 10Li3(l -z) + 8Li3(z) 
3 2 5 

- 8Li2(z) log z + 10Li2(z) log(l -z)- 4C2 log z - IOC2 log(l - z) - log3 Z 
+ 51og2(l-z)logz-8C3 + 0(e')l i, 
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Yf(z,e) = (l-z)- 



-l-3e 



1 



^31ogz + 



1 



Li2 (z) + log^ z + C2 + 6 log 2; 



(A.IO) 



+ Li3(l -z) + 2Li3(z) + U2{z) log(l - 2;) - C2 log(l -z) + 8C2 log z 

+ - log2;log2(l -z)- 2(3 + 2Li2(z) - 21og2 z - 2(2 + e[6C2Li2(z) - 2Li2(z) 



- 10Li4 1 - 



9Li4(l -z)- 6Lu{z) - 2U2{z) log^ z - -U2{z) log^{l - z) 



+ 8Li3(z) log z - Li3(l - z) log(l -z)- 10Li3(z) log(l - z) - 6(2 log' z 
+ -C2 log'(l -z) + IOC2 log(l - z) log z + 8C3 log z + IOC3 log(l - z) 



7 
12 



log z 



1 



+ - log(l - z) log^ Z-- log^(l - z) log z - - log^(l - z) log' z - -Q - 2Li3(l - z) 
- 4Li3(z) - 2Li2(z) log(l -z) + 2(2 log(l -z)- I6C2 log z-\ogz \og'^{l - z) + 4(s 



Yf(z,6) = (l-z)- 



l-4e 



^21og2:+ - 



2Li2(z)-2C2 + 41ogz +2Li3(l-^) (A.ll) 
2 



+ 4Li3(2;) + 2Li2(z) log(l - z) - 2(2 log(l - 2;) + - log'^ z + log^(l - z) log 2; 

o 

- 4C3 + 4C2 - 4Li2(z) + e[l2C2Li2(z) - 8U2{zf - 20Lu (1 - 



- 22Li4(l - z)- 20Li4(z) - 8Li2(z) log' z + Li2(;z) log'(l - z) + 24Li3(z) logz 
+ 2Li3(l - z) log(l -z)- 20Li3(z) log(l -z)- 2(2 log' z - (2 log'(l - z) 

+ 2OC2 log(l - z) log z + 8(3 log z + 2OC3 log(l - 2) - ^ log^ ^ 



10 2 
+ — log(l - z) log^ z + - log^(l - log 2; - 5 log'(l - z) log' z + 4(2 
3 3 



2 

4Li3(l -z)- 8Li3(z) - 4Li2(z) log(l - z) + 4(2 log(l - ^) - ^ log^ z 



21og'(l-z)logz + 8C3 +0(e' 



A. 3 Results for the double-real master integrals 



Xf(z,e) = zlogz- J(l-z)(z'-5z-2) +e[^(l-z)(z'-5z-2)log(l-z) (A.12) 
+ 42;Li2(2;) - 4z(2 - ^2; log' z - ^z{z + 2) log 2; - ^(1 - z) {4z^ - 17 z - 5) 



+ e'| -8z\ogz- z){z^ -hz-2) log'(l 
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+ 16zC2-16zLi2(z)- (l-z2)logz + -(l-z)(4z2-l7z-5) log(l - z) 

1 •J 

- 16zLi3(l - z)- 8zUs{z) + 2zU2{z) log z + SzCs + 2zC2 log ^ + 772; log"^ z 

o 

- + 4^ + 1) Li2(z) - ^(2-2 + 1) (^^ - 8z + 1) C2 + ^z{z + 2) log2 z 

- ^z{5z + 8) log z - ^(1 - z) {U8z^ - 1625z - 281) } 

+ [yzlogz + ^(1 - z) {z^ - 5z - 2) ] log3(l - + [32zLi2(z) - 32^2 

+ 4z log^ z + 4(2z + 1) log z - ^(1 - z) (4z^ - 17z - 5) 1 log^(l - z) 

9 J 

+ 64zUs{l - z) + 32zLi3(z) - 8zLi2(z) log z - 32zC3 - 24zC2 log z 

16 1 

- —z Iog3 z + 4{z'^ + 4z + l) U2{z) + 7: {I- z^) log2 z 

+ ^{2z + 1) (z^ _ 8z + 1) C2 - ^ (1 - z') log z 

+ ^(1 - z) (448^2 - 1625z - 281) log(l - z) - 16zC2Li2(z) 



+ 32zLi4 1 



- 32zLi4(l - z)- 28zhu{z) - zLi2(z) log^ z 

- 8zLi3(l - z) log z + 2zLi3(z) log z + -^zC,l + 15^(2 log^ z + lOzCs log z 

o 

+ log^ z + 4 (z^ + 4z + 1) Li3(l - z) + (5^2 + 8z - 1) Li3(z) 

- (2^2 + 2^-1) Li2(z) log z + ^ (-8^2 + 33^2 - ASz - 13) Ca - z{z + 2)C2 logz 

- + 2) log=^ z + ^ (-5^2 - 16z - 5) Li2(z) 

+ (-40z^ + 255^2 + 2Az - 5) C2 + lz{5z + 8) log^ z - lz{3z + 4) log z 
18 8 8 



-, 1^77(1 - z) (I0496z^ - 33385Z - 1897) ) + O(e^) , 
1296 J 



Xf (z, e) = - logzlog(l -z)- U2{z) + + C2 + ^ - 1 + e{41ogzlog2(l - z 



(A.13) 



+ 



4Li2(z) - - log2 z - 4C2 + 2 log z + 4(1 - z) log(l - z) + 6Li3(l - z) 



+ 7Li3(z) - 2Li2(2;) log z - 8(2 log z - - log^ - 7(3 + 2Li2(2;) - log^ z 

- 2(2 - 4(1 - z)} + - 81ogzlog3(l - z) + [ - 8Li2(z) - Slog^ z + 8(2 

- 8 log z - 8(1 - z)l log2(l -z)+\- 24Li3(l -z)- 28Li3(z) + 2Li2(z) log z 
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+ 3OC2 \ogz + — log^ z + 28C3 - 8Li2(z) + log^ z + 8C2 + 16(1 - z] 



log(l - z) 



- 19 - 



+ 10C2Li2(z) - SU2{zf - 22Li4 (^1 - - J + 6Li4(l - z) + 5Lu{z) + U2{z) log^ z 

+ 2Li3(z) log z - y C2 log' z - 7C3 log ^ - ^ log^ ^ - 9Cl " 12Li3(l - z) 

2 o 

- 14Li3(2;) + 4Li2(z) log z + IOC2 log 2; + - log'^ z + 14(3 + 4(1 - z)C2 



16(l-z)}+0(e3 



X|^(2, e) = ^ [ - logz + z - ij + - [4(1 - z) log(l -z)- 4Li2(z) + log^ z + 4(2 (A.14) 



+ 4 log z 



+ 



log z- 8(1 -z) log2(l-z)+ 16Li2(z) - I6C2 log(l-z) 



+ 16Li3(l -z) + UUsiz) - 2Li2(z) log z - 6C2 log z - - log^ z - 12(3 
+ 16Li2 (z) - 4{z + 3)C2 - 4 log^ z - 4 log z - 4(1 - z) 



64 32 

log-z + ^(l--2) log^(l-z)+ -32Li2(z)-81og^z + 32C2 



3 ' 3 

2/ 



32 log z log^(l-z)+ - 64Li3(l -z) -48Li3(z) + 8Li2(z)logz + 40C2logz 



+ 81og3 z + 48C3 + 16(2 + 3)C2 - 64Li2(z) + 16(1 



log(l -z) + 16C2Li2(z) 



- 48Li4 ^1 - + 16Li4(l - z) + 26^4(2) + ^2(2) log^ z + SLi^il - z) \ogz 

5 132 

- I8C2 log' z - IOC3 log z - - V z - — Cl - 64Li3(l -z)- 48Li3(z) 
+ 8Li2(z) log z - 16(z - 4)C3 + 24C2 logz + - log=^ z - ISLhiz) + 41og' z 
+ 16C2-16(l-z)}+0(e'), 



Xf(z,e) = 



{l-z) 



-4e 



+ 



(1 - z)e^ 



2 - 4e + (2 - 8C2)e' + (I6C2 - 32C3)e^ + 0{e^) 
- 2 log z- 2(1 - z 

+ 4C2 - 4Li2(z) + 81ogz + 4(1 - z)} + - 16(1 - z) log'(l - z) 

16Li2(z) - log2 z - I6C2 - 16 log z - 16(1 - z)] log(l - z) + 10Li3(l - z) 



(A.15) 



+ ^^-4^{[41ogz + 8(l-z) log(l-2) 



+ 



+ 6Li3(2;) - 4Li2(z) log z + 6(2 log z + 



\og^ z 



-6C3 + 16Li2(z)-8(2 + l)C2 



32 64 

_logz+_(l-2)j log3(l-z) 



-81ogz-2(l-z)} + ^{ 

+ [ - 32Li2(z) - log' z + 32C2 + 32(1 - z)] log'(l - + [ - 40Li3(l - z) 
- 24Li3(z) + 6Li2(z) log z + 2C2 log z - log^ z + 24(3 - 64Li2(z) + 32(z + 1)C2 
+ 41og' z + 161og z + 8(1 - z)] log(l -z) + 26C2Li2(z) - 5Li2(z)' 
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+ 6Li4 1^1-- j + 30Li4(l - z) + 26Li4(z) + 3Li2(z) log^ z - 10Li3(l - z) log z 

- ULhiz) log z - 2C2 log2 z + 2OC3 log z + ^ ^Cl - 40Li3(l - z) 

- 24Li3(z) + 16Li2(z) log z - 8{4z - 7)C3 - 24C2 log 2 - ^ log^ z + 16^(2 



- 16Li2(z)} +0(€ 



Xf(^,e) 



(1-z) 



-4e r 



1 + 3e + (4C2 - 4)6^ + (-I2C2 + I6C3 + 2)6^ 
+ (yC| + 16C2-48C3)e^ + 0(e^) 
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1 - z 



+ 



4(1 - zf log(l - z) + log ^ - (1 - z){3z - 4) 



(1 - z)ze2 

+ - 8(1 - ^? V(l - ^) + 4(1 - zK3z - 4) log(l -z)+ 4Li2(z) 

+ 4(z - 2)<2 - 4 log z - 4(1 -zf] + ^^-^-^{^(1 - zf log3(l - z) 

8 log z - 8(1 - z){3z - 4)1 log2(l -z)+\- 16Li2(z) - 16(z - 2X2 
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+ 16(1 - zf\ log(l -z)- 10Li3(l -z) + 4Li3(z) - 2Li2(z) log z 

+ 4{2z - 3){2z - 1)C3 - 6C2 log z - ^ log3 z - 16Li2(^) - 4^(3^ - 7)C2 + 4 log z 

- 2(1 - z)z] + 7Y^{ - f (1 - log^(l - ^) + [f log ^ 
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+ y (1 - z){3z - 4)J log'^(l -z)+ ^32Li2(z) + 32(z - 2)zC2 - 91og^ z + 321ogz 

- 32(1 - zf log2(l -z)+ 40Li3(l - z) - 16Li3(z) - 10Li2(2) log z 

- 16(2z - 3)(2z - 1)C3 + 26C2 logz + ^ log=^ z + 64Li2(z) + 16z{3z - 7)C2 



+ 8(1 - z)z log(l -z) + 2C2Li2(z) - 9U2izf - 2Lu 



26Li4(l -z) 



- 26Li4(2:) + 4Li2(z) log^ z - 18Us{l - z) log z - 2lAz{z) log z 
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+ 21og^z + 81ogz + 8(l - z)} + - 16 log 16(1 -z) \og^{l-z) 



+ 
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+ 4Li2(-z) log z + I2C2 log z - 44C2 log(l + z) + ^ log^ z-^ log^(l + 
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+ ^ log^ 2 log 2 - ^ log^ (1 + z) log2 2-2 log2 2 log(l + 2) log ^ + ^Cl 
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+ 16C2log2 



16^2 2, 



Xf (2, 6) 



(2 + l)e2 



2 log 2 + 2^ - 1 
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+ 4Li2(-2) - 8Li2(2) - log^ 2 + 41og(l + 2) log 2 + IOC2 + 8 log 2 

+ [l6l0g2 - 8 (1 - 2^) ] log2(l - 2) + [ - 16Li2(-2) + 32Li2(2 

+ 8 log2(l + 2) - 16 log 2 log(l + 2) - 16 log 2 log(l + 2) - 4OC2 



+ 81og2 2 
-8Li3 



log(l - 2) + 16Li3 



1-2 



+ 16Li3(l - 2) - 12Li3(-2) 



1 \ /l-2\ /1 + 2 

+ 16Li3 —— + 16Li3 



1 + 2 



1 + 2 



BLi2(— 2) log 2 
4. 



(A.19) 



- 4Li2(2) log 2 + IOC2 log 2 - 2OC2 log(l + 2) + - log^ 2 - - log^(l + 2) 

o o 

16 

+ 2 log(l + 2) log2 z + 8 log^ 2 log(l + 2) - I8C3 + I6C2 log 2 - — log3 2 

- 16Li2(-2) + 32Li2(2) - 4 (2^ + 9) C2 + 41og22 - 16 log(l + 2) log 2 - 81og2 

- log 2 + I (1 - 2^) ] log3(l -z)+ [32Li2(-2) 



2 



- 64Li2(2) + 8 log^ 2-32 log^(l + 2) + 32 log(l + 2) log 2 + 64 log 2 log(l + 2) 



+ 8OC2 - 321og^2 - 64 log 2 

1 



log2(l -z) + 



64Lia 



1 - 2 



64Li3(l - 2) 



+ 48Li3(-2) + 32Li3 



1 + 2 



64Li3 



1 - 2 

1 + 2 



64Li3 



1 + 2 



-23- 



- 8Li2(-z) log z + 16Li2(z) log z + 4C2 log z + 8OC2 log(l + z) - 8 log^ z 
+ — log^(l + z) + 161og(l + z) log^2; + 41og^(l + z) log2; + 41og^21ogz 

- 32 log^ 2 log(l + z) - 8 log 2 log(l + z) log z + 72(3 - 64(2 log 2 + — log^ 2 

3 

+ 64Li2(-z) - 128Li2(z) + 16 [z^ + 9) C2 - 321og2(l + z) + 641ogzlog(l + z) 
+ 64 log 2 log(l + z) - 32 log^ 2 + 16 (1 - z^) log(l - z) + 12Li4 (^1 - ^ 
+ 4Li4 (1 - z^) - 38C2Li2(-z) + 32C2Li2(z) + 2\A2{,-zf + 24Li2(-z)Li2(z) 



- 48Li4 ^1 - ^ j + 48Li4(l - 2) - 34Li4(-z) + 8Li4(z) + 24Li4 (^j^ 

+ 48Li4 (ytI) ~ ^^^'^ (^Tt) ^ '^^^'^ (^) ~ ^^''^"''^ ^ 
+ 2Li2(2;) log^ z + 4Li2(-2;) log(l + z) log z + SLig (-y^) log ^ 

+ 8Li3(l - z) log 2 + 16Li3(-2) log z - 12Li3(2) log z + 12Li3 ( log , 



+ 8Li3 J ^ ^ ^^'^ J ^ ^ 48Li3(z) log(l + z) - 13C2 log' z 

- 24C2 log2(l + z) - 28C2 log(l + 2) log z - 44C3 log z + 8C2 log 2 log z 

21 56 
+ 36C3 log(l + z) + — log^ 2; + 41og^(l + 2;) - Y log(l + z) log^ z 

- ^ log^(l + z) log z - ^ log^ 2 log z + 20 log2(l + z) log^ z 

+ 41og2 2 log(l + z) log z - ^Cl - 64Li3 (^^) - 64Li3(l - z) + 48Li3(-z) 
+ 32Li3 \-^—- \ - 64Li3 ( ) - 64Li3 ( ^^"^ ) - 32Li2(-2) log z 



\ + z) \\ + z) \ 2 

X6 

+ 16Li2(z) log z - 8 (2^2 - 11) C3 - 4OC2 logz + 8OC2 log(l + ^) - — log^ ^ 

+ ^ log2(l + z) - 8 log(l + z) log' z - 32 log' 2 log(l + - 64(2 log 2 + ^ log^ 2 

+ 16Li2(-z) - 32Li2(z) - 41og2 z + 161og(l + z) log z + 40(2 - 16 (l - z^) } 



log z + z — 1 



1 

+ - 
e 



4(1 - z)log(l 



4Li2(z) + log'z + 4C2 (A.20) 



+ 41ogz + 81ogz-8(l-z) log'(l-z)+ 16Li2(z) - I6C2 log(l - z) 

+ 16Li3(l - z) + 12Li3(z) - 2Li2(2;) log z - 6(2 log -2 - ^ log^ z - I2C3 

+ 16Li2(z) - 4(z + 3)C2 - 4 log' z - 41ogz - 4(1 - z) 



- 24 - 



+ e{ [ - y logz + y (1 - z)\ log3(l -z)+[- 32U2{z) - 81og2 Z + 32C2 

- 32 log z log^ {l-z)+ - 64Li3 {1 - z) - ASU^ (z) + 8U2 {z)\ogz 

+ 4OC2 log z + 8 log3 z + 48C3 + 16(z + 3)C2 - 64Li2(z) + 16(1 - z)] log(l - z) 



+ 16C2Li2(z)-48Li4 1 



+ 16Li4(l -z) + 26Li4(z) + Li2(z) log^ z 



+ 8Li3(l - z) log z - I8C2 log^ z - lOCa log z - - log^ z 



132 



Cl - 64Li3(l - z) 



48Li3(z) + 8Li2(z) log z - 16{z - 4)C3 + 24(2 log z + - log^ z - 16Li2(z) 

o 



+ 41og2 z + I6C2 - 16(1 - 2)} + 0(e2) , 



Xfo(^,e) 



Z€^ 



log ^ — (1 — z)z 



+ 



Z€ 



4(1 - z)z log(l - z)- 4Li2(-z) 



(A.21) 



4Li2 {z) + 2 log^ z - 4 log(l + z) log 2; + 2C2 + 4 log z 



+ 



Slogz - 8{1 - z)z log2(l-z)+ 16Li2(-z) + 16Li2(z)-81og2(l + z) 



+ 16 log z log(l + z) + 16 log 2 log(l + z) - 8C2 - 8 log^ 2 log(l - z) 



- 16Li3 

- 16Li3 



1 - z 
2 

1 



+ 32Li3(l - z) + 14Li3(-z) + 20Li3(z) + 12Li3 

1 + z 



l + z 



16Li3 



6Li2(-2:) log z - 2Li2(-z) log z — IOC2 log z 



+ 22C2 log(l + z)-2 log^ z + - log=^(l + z) + log(l + z) log2 z-8 log^ 2 log(l + z) 
16 

- 4C3 - I6C2 log 2 + — log3 2 + 16Li2(-z) + 16Li2(z) - 4 {z'^ - z + 2) (2 

- 8 log^ z + 16 log(l + log 2; - 41og2;- 4(1 - z)z^ + -| 

32 1 r 

+ — (1 - z)z log3(l -z)+ - 32Li2(-z) - 32Li2(2) - 161og2 z + 321og2(l + z) 

- 32 log(l + z) log z - 64 log 2 log(l + z) + I6C2 + 32 log^ 2-32 log z] log^(l - z) 



64 

- Y log 



+ 64Li3 
+ 64Li3 



1 - z 
2 

1 - z 
l + z 



128Li3(l - z) - 56Li3(-2) - 80Li3(z) - 48Li3 
l + z 



l + z 



+ 64Li3 



68 



+ 8Li2(-z) log z + 8U2{z) log z + 48(2 log z 



- 88C2 log(l + z) + — log-" Z-- log''(l + z)-20 log(l + z) log^ z 

- 4 log^(l + 2;) log z - 4 log^ 2 log z + 32 log^ 2 log(l + z) + 8 log 2 log(l + z) log z 
+ I6C3 + 64C2 log 2 - ^ log=^ 2 - 64Li2(-z) - 64Li2(z) + 16 (z^ - z + 2) C2 

+ 32 log^ (1 + z)-64 log 2 log(l + z) - 64 log 2 log(l + z) + 32 log^ 2 



-25- 



+ 16(1 - z)z^ log(l -z)- 14Li4 (^1-^^- 2Li4 (l - z^) + 34C2Li2(-z) 
+ 16C2Li2(z) + 2Li2(-z)2 - 16Li2(-z)Li2(z) - 2AU4 (^1 - + 56Li4(l - 

+ 14Li4(-z) + 70Li4(z) - 2U4 (^YTz) ~ ^^^^^ (^T^) ^^^^^ ( 

62Li4 ( + 2Li2(--z) log^ z + Li2(z) log^ 2; + 4Li2(-2;) log(l + z) log z 
log z + 16Li3(l — z) log z + 8Li3(— z) log z + 4Li3(2;) log z 



z) 

z-V 



- 62Li4 




(z + 1 


- 8Li3 


^) 




+ lOLia 







log z - 8Li3 ^ ~ ^^^^ (~^) ^ 

- 32Li3(z) log(l + z) - 24C2 log2 z + I6C2 log^ll + + 32(2 log(l + z) log z 
+ 2C3 log z - 8C2 log 2 log z - 24C3 log(l + z) - 9 log^ z - ^ log^(l + z) 

+ — log(l + z) \o<i z + 10 log-'^(l + z) log 2; + - log'^ 2 log z - — log2(l + z) \o<i z 



- 41og^ 21og(l + z) \ogz- -CI + 64Li3 ) ' 128Li3(l - z) - 56Us{-z) 

- 80Li3(z) - 48Li3 ( + ^^^'3 (tT^ ) ^ ^'^^'^ (^^) ^ ^ 
+ 8Li2(z) log z-16{z^ - z-l)C3 + 4OC2 log z - 88C2 log(l + z) + 8 log^ z 

- ^ log3(l + z) - 4 log(l + z) log2 z + 32 log^ 2 log(l + z) + 64(2 log 2 - ^ log^ 2 

- 16Li2(-z) - 16Li2(z) + 81og2 z - 161og(l + z)\ogz + 8(2 - 16(1 - z)z^ 
+ 0(e'), 



Xfl„(^,6) 



- ^ (2 - 3)^2 log z - -^(1 - z) (5^2 - 22z + 5) 
+ [ - 1(1 - z)3 logz + ^(1 - z) (5^2 - 22z + 5) ] log(l - z) 



(A.22) 



^{z + 1) (z^ - 4z + 1) Li2(z) + ^(^ + 1) (-2' - 4z + 1) C2 + ^{z - 3)^2 V z 

—z ihz^ 
12 ^ 



7z + 4) logz - ^(1 - (205z^ - 6622 + 205) | 
+ e^l [ - \{Zz - l)logz - ^(1 - z) (5^2 - 22z + 5) 1 log2(l - z) 



+ 



|(z + 1) (z2 -42 + 1) Li2(2) - ^(^ + 1) {z^ - 4z + 1) C2 + ^(1 - zf \oi z 



1 

6' 



- ^(1 - z) (5^2 - 6z + 5) logz + -^(1 - z) (205^2 - 662z + 205) log(l - z) 
+ ^(2 + 1) [z^ - 4z + 1) Li3(l - 2) + ^ (52^ - 15^2 + 3z - 1) Li3(2) 



-26- 



+ i (-2^3 + -3z + l) U2{z) log z-^{z- 3)z2C2 log z-j^{z- 3)z^ log3 Z 

+ ^ {-5z^ + 15^2 - 3z + 1) Cs - ^(-z + 1) (5^^ -8z + 5) U2{z) 

+ ^z (5^2 -7z + 4) log2 z + (5^3 + 45^2 - 63z + 25) C2 

- -^z (27z^ - 25z + 28) logz - -^(1 - z) (6365^^ - 17158z + 6365) ) 
24 1296 J 

8 81 
- - (z^ - 3z^ -9^ + 3) logz + —(1 - z) {bz"^ - 22z + 5) log3(l - z) 
y ^1 



+ 6^ 



+ 



-(z + 1) (z^ - 4z + 1) Li2(z) + -(z + 1) - 4z + 1) C2 



1 2 
+ - (-3^3 + dz"^ -3z + l) log^ z + - {Az^ -7z + 5) log z 

- -[-(1 - z) (205z^ - 662^ + 205) ] log2(l - z) 



+ 



16 2 

— + 1) {z^ -Az + l) Li3(l - z) + -(z + 1) (z^ -4z + l) U2{z) log , 



- ^ (5/ - 15^2 + 3z-l) Li3(z) + ^ (5^3 - 15^2 + 3z - 1) Cs 

+ - (7^3 - 21^2 + 3z-l) C2 log z + ^ {9z^ - 27z^ - 2\z + 7) log^ z 
3 18 

+ h^z + 1) (5^2 - 8z + 5) Li2(z) + ^(1 - ^) (5^^ -6^ + 5) log^ z 

- \ (5^3 + 45^2 - 63^ + 25) C,2-\-{\-z) {21 z^ - 26z + 27) log z 
9 12 

+ ^(1 - z) (6365^2 - 17158Z + 6365) log(l - z) 

- ^{z + 1) (z^ - 4z + 1) Li4 (^1-^^ +^{z + l) (z^ - 4z + 1) Li4(l - z) 

+ l{z + l) {z^ - 4z + 1) hu{z) + \ {bz^ - 15^2 + 3^-1) C2Li2(2) 

+ ^ (2^3 - 6^2 + 3z - 1) Li2(z) log2 z-\{z'^- 3z^ + 6z - 2) Li3(l - z) log z 

+ \ {-2z^ + 6^2 - 3z + 1) Lisiz) log z + (1 - z)3Li2(z)2 

+ — (-49^3 + 147^2 + 57^ - 19) C| + J {-7z^ + 21^^ ^ 24z - 8) C2 log^ z 

15 D 

+ ^ {3z^ - 9z^ + 24z - 8) Ca log z + ^ (-15z''^ + 45z^ + 48z - 16) log"^ z 

2 1 
+ -{z + 1) (5^2 - 8z + 5) Li3(l -z) + - {25z^ - 39^^ + 27z - 5) Li3(^) 

o D 

+ ^ {-Wz^ + ISz^ - 152; + 5) Li2(2) log z-^z (5^2 -Iz + A) (2 log 2 

- ^z (5^2 -7^ + 4) log3 z + ^ (-115^3 + 333^^2 _ 297^ + 55) ^3 
72 18 

-^{z + 1) (9^2 - 8z + 9) Li2(z) + ^z {27 z^ - 25z + 28) log^ z 



-27- 



+ (-167z^ + 1761^2 - 1707z + 653) C2 - (153^^ - 131z + 164) logz 
108 48 



7776 



(l-z) (35537^2 - 86302Z + 35537) } + 0(e^ 



Xf2(z,e) = -(1 - z)logzlog(l -z)-{l- z)Li2{z) + (1 - z)C2 + ^(1 - z)log'z (A.23) 
+ 21ogz - ^(1 - z) (z^ -2z- 23) + e|4(l - z) logzlog2(l - z) 



+ 



4(1 - z)U2{z) - 4(1 - z)C2 + -iz-l) log2 z + 2(1 - z) log z 



+ -(1 - z) (z^ -2z- 23) log(l - z) + 6(1 - z)Li3(l - z) + 7(1 - z)U3{z) 

- 2(1 - z)Li2(z) log z - 7(1 - z)C3 - 5(1 - z%2 log z + ^(z - 1) log3 z 

- 2{z - 5)Li2(z) + 2{z - 5)C2 + (-^ - 2) log^ z - 2zlogz 

- ^(1 - z) (5^2 - lOz - 103) } + - 8(1 - logzlog^(l - z) 
+ [ - 8(1 - z)U2{z) + 8(1 - z)C2 - 5(1 - z) log2 z + 8{z- 3) log z 

- |(1 - z) (z^ - 2z - 23) j log2(l -z)+\- 24(1 - z)Li3(l - z) 



28(1 - z)U3{z) + 2(1 - z)U2iz) log z + 28(1 - z)C3 + 30(1 - z)C2 log z 

8(z-5)C2 + (l-2)log^2 



25 q 

+ — (1 - z) log3 z + 8(z - 5)Li2(z 



- 4(1 - z) log z + - (1 - z) (5^2 - lOz - 103) J log(l - z) 



+ 10(1 - z)C2U2iz) - 3(1 - z)U2{zy - 22(1 - z)U4 y--- 

+ 6(1 - z)Li4(l -z) + 5(1 - z)U4iz) + (1 - 2)Li2(z) log^ z 

15 

+ 2(1 - z)Li3(z) log z - 9(1 - 2)C| - y(1 - ^)C2 log' z - 7(1 - 2)C3 log z 
19 

-—(l-z) log^ z + 4(3z - ll)Li3(l -z) + 2{7z - l5)Lh{z) 

- 4{z - 2)Li2(z) log z - 2{7z - 15)C3 - 2(5z - 7)C2 logz + ^(3 - 2z) log^ z 

- 4(2; + l)Li2(2) + \ {-z^ + Sz^ + 332 - 11) C2 + -zlog' z - 6z\ogz 



— (1 - z) {Uz^ - 28z - 229) } + 0{e^) , 



Xg{z, e) = -2Li3(z) + Li2(z) log z + C2 log z + J log" z + 2(3 + 1(1- zf (A.24) 

6 2 

+ e{ log2 z log2(l -z)+ [8Li3(z) - 2Li2 {z) log z - 6C2 log z - ^ log^ z - 8(3 



- 2(1 - zf log(l -z)- 2C2Li2(z) + Li2(z)2 + 8Li4 



04) 



-28- 



+ 8Li4(l - z) + 8Lu{z) + -U2{z) log^ z - 2Li3(l - z) log z - 5Li3(z) log z 

+ ^C2 log^ z - 3C3 log z + ^ log^ - y Cl + 8Li3(^) - 4Li2(z) log z - 4(2 log z 



- Vz-8C3 + 2(1-^)2} +0(e2), 



Xf4(^,6) 



log^ z 1 



(1 - zf + 2(1 - zf log(l - z) + 4Li3(z) - 4C2 log z (A.25) 



--log3z-4C3-21og^z-(l-z)2 + e| -41og2z - 4(1 - log2(l - z 
+ [ - 16Li3(z) + I6C2 log z + ^ log^ z + I6C3 + 4(1 - zf^ log(l - z) 

- 16Li4 ^1 - i j - 16Li4(l -z)- 2Li4(z) + 2Li3(z) log z - 3(2 log^ 

- ^ log^ z + ^Cl - 16Li3(z) + I6C2 log z + 2 log3 z + I6C3 + 2(1 - zfC2 
+ 21og2z-4(l-z)2}+C»(e2), 



{z + l)en 



logz + -{z^-l) 



+ 



(z + l)e 



2 (1 - log(l - z) 



(A.26) 



+ 4Li2(-z) - 4Li2(z) - log^ z + 41og(l + z) log ^ + 6(2 + 4 log z 



+ 



^1 81ogz- 4(1-^2) log2(l-z)+ - 16Li2(-z) + 16Li2(z 



2; + 

+ 81og2(l + 2;) - 161ogzlog(l + z) - 161og21og(l + z) - 24C2 



+ 81og^2 



+ 16Li3 



log(l -z) + 16Li3 



1 - z 



8Li3(-z) - 8Li3(z) + 16Li3 



1 - z 



1 + z 



1 + z 



+ 8Li2(-2;) log z + I2C2 log z - I6C2 log(l + z) + log^ z 



- log^(l + z) + 4 log(l + z) log2 z + 8 log2 2 log(l + z)- 14(3 + I6C2 log 2 
^^^y ^ - 16Li2(-z) + 16Li2(z) - 2 (^2 + 11) C2 + 41og2 z - 161og(l + z) logz 



41ogz-2(l-z2)} + _A_| 



64 , 16 2 

— log z -\ (1 — z 

3 ^ 3 ^ 



log3(l - z) 



+ 



32Li2(-z) - 32Li2(z) + lOlog^ z - 321og^(l + z) + 321og(l + z) log z 



+ 64 log 2 log(l + z) + 48C2 - 32 log^ 2-32 log z log^ (1 - z) 
+ [ - 64Li3 (^y^) + 32Li3(-z) + 32Li3(z) - 64Li3 (^^^ 
- 64Li3 ( j + 4Li2 (z) log z - 2OC2 log z + 64(2 log(l + z) - y log^ z 



-29- 



+ — log3(l + z) + 16 log(l + z) log2 z - 32 log2 2 log(l + z) + 66(3 - 64(2 log 2 

o 

64 

+ — log=^ 2 + 64Li2(-z) - 64Li2(z) + 8 (z^ + ll) (2 - 321og2(l + z) 
+ 64 log z log(l + z) + 64 log 2 log(l + z) - 32 log^ 2 + 8 (l - z^) log(l - z) 
+ 8U4 (^1-^^+ 8Lu (1 - z'^) - 48C2Li2(-z) + 12C2Li2(z) + 2U2{zf 
+ 32Li2(-z)Li2(-2) - 24Li4(-z) - 10Li4(z) + 64Li4 f + 64Li4 

6AU4 I I + 64Li4 ^ 



z + 1 



z + l 



1 + zJ \l + z 

12Li2(-z) log2 z + 2Li2(z) log2 z 



+ 4Li3(l - z) log z + 24Li3(-2;) log z - 14Li3(2;) log z + S2U3 



1 



l + z 



logz 



+ 64Li3(z) log(l + z) + 5C2 log^ z - 32C2 log2(l + z)- 32(2 log(l + z) log z 

- 24C3 log z + 48C3 iog(i + + Y log^ ^ + y iog^(i + ~ y ^^g(i + ^) ^°g^ ^ 



16 log^(l + z) log z + 16 log^ (1 + z) log^ z 



+ 32Li3(-z) + 32Li3(z) - 64Li3 



1 - z 
l + z 



64Li3 ( -32Li2(-z)logz 

32 



- 8 (z^ - 8) C3 - 48C2 log z + 64C2 log(l + z) - 4 log^ z + — \og^{l + z) 

o 

64 

- 16 log(l + z) log^ z - 32 log^ 2 log(l + z) - 64C2 log 2 + — log=^ 2 + 16Li2(-z) 

3 

- 16Li2(z) - 4 log2 z + 16 log(l + z) log z + 24(2 - 8 (l - z^) } + 0{e^) , 



Xf6(^,6) = 



1 - Z 



+ 



2ze' z{z + l)e2 



- 2 (1 - z^) log(l - z) + z\ogz-2{l- z^) 



(A.27) 



+ 



z{z + l)e ~ - ^) + 8 (1 - z") log(l -z) + 2(1 - z)Li2(-z) 

+ 4zLi2(z) + (2z^ - 5z - 1) C2 - log^ z + 2(1 - z) log(l + z) log z - 4z log z 



8z log z 



16 (1 - z' 



z(z + 1)1 3 
Iog2(l - z) + [ - 8(1 - z)Li2(-z) - 16zLi2(z) 

- 4 (2z2 - 52 - 1) (2 + (z + 1) log2 z + 4(1 - z) log2(l + z) + 4(1 - z) log^ 2 

- 8(1 - z) log(l + z) log z - 8(1 - z) log 2 log(l + z) - 8 (l - z^) ] log(l - z) 

1 - z' 



+ 8(1 - z)Li3 



10(z + l)Li3(l - z) - 4(z + 3)Li3(-z) 



- 2(3z + 7)Li3(z) + 8(1 - z)Li3 ( ^—- \ + 8(1 - z)Li3 ^ ^ ^ ^ 



1 + zy ■ V 2 

+ 8Li2(-z) log z + 4(z + l)Li2(z) log z + {Sz^ + llz - ll) (3 



-30- 



- A{z - 2)C2 log z - 8(1 - z)C2 log(l + z) + 8(1 - z)C2 log 2 + ^(z + 4) log^ z 

- ^(1 - z) log3(l + z)-^il-z) log3 2 + 2(1 - z) log(l + z) log2 z 
+ 4(1 - z) log2 2 log(l + z) - 8(1 - z)U2{-z) - 16zLi2(z) 

- 4 (22;^ - 52; - 1) C2 + 4 log^ z - 8(1 - z) log(l + 2;) log z + 4z log 

+ Y (1 - log'(l - ^) + [y^logz + ^ (1 - 2^) ] V(l - 2) 



+ 



16(1 - z)U2{-z) + 32zLi2(z) + 8 {2z^ - 5z - l) C2 + (^ + H) log 
- 16(1 - z) log^(l + z)- 16(1 - z) log^ 2 + 16(1 - log(l + z) log z 
+ 32(1 - z) log 2 log(l + + 32z log z + 16 (1 - z'^) 1 log2(l - z) 



'z 
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32(1 - z)Li3 ( i-^ ) + 40(z + l)Li3(l -z) + 16{z + 3)Li3(-2) 
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+ 32(1 - z)C2 log(l + 2) - 32(1 - z)C2 log 2 + ;^(-13z - 45) log^ 2 

16 32 
+ —(1 - z) log^(l + z) + — (1 - 2) log=^ 2 + 8{z + 3) log(l + z) log2 z 

o o 

- 16(1 - 2) log2 21og(l + 2) + 32(1 - 2)Li2(-2) + 64202(2) 

+ 16 (22^ - 52 - 1) C2 - 4(2 + 1) log^ 2 - 16(1 - 2) log2(l + 2) - 16(1 - 2) log^ 2 
+ 32(1 - 2;) log(l + 2) log z + 32(1 - z) log 2 log(l + 2) log(l - z) 
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- 2(132 + 7)C2Li2(2) - 8(2 + l)Li2(-2)^ + (52 + 7)Li2(2)^ + 32Li2(-2)Li2(2 

- 10{z + l)Li4 ( 1 - - ) - 54(2 + l)Li4(l -z) + 4{7z + l)Li4(-z) 



2(92 + 14)Li4(2) + 64Li4 ( 7^ ) + 32(1 - 2)Li4 — - 
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- 32(1 - z)Li4 (^T^J + 64Li4 + 2(5^ - l)Li2(--2) log' z 

- (32 + l)Li2(2) log2 2 - 16(2 + l)Li2(-2) log(l + z) log 2 
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+ 5C2 log^ z - 32C2 log2(l + z)- 16(z + 3)C2 log(l + 2) log 2 + 6(32 - 1)C3 log z 

1 16 
+ 48C3 log(l + 2) + —{2lz + 79) log^ 2 + y log^(l + 2) 



- 31 - 



2 16 
--{5z + 27) log(l + z) log3 z + —{z-2) log^ (1 + z) log z 



+ 8(1 - z) log2(l + z) log2 2 - 32(1 - z)U3 



l-z 



+ 40(z + l)Li3(l-z) 



+ 16(z + 3)Li3(-z) + 8{Sz + 7)Li3(z) - 32(1 - z)U3 



l-z 
l + z 
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- 32Li2(-z) log z - 16{z + l)Li2(z) log z 
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4 o 16 o 

- 32(1 - z)C2 log2 - - (z + 4) log3 z + —{l-z) log3(l + z) 

32 

+ —{l-z) log^ 2 - 8(1 - z) log(l + z) log^ z - 16(1 - z) log^ 2 log(l + z) 
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Xf7(z,e) 



(l-z) 



-4e 



2 - 4e + (2 - 8C2)e^ + (I6C2 - 32(3)6^ + 0(e^ 
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(1 - z)e2 L 



-21ogz-2(l -z) 



+ (fZ^ { [4 log z + 8(1 - z)] log(l - z) 

+ 4C2 - 4Li2(z) + Slog z + 4(1 - z)} + - 16(1 - z) log2(l - z) 

16Li2(z) - log2 z - I6C2 - 16 log z - 16(1 - z)l log(l - z) + 10Li3(l - z) 
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+ 6Li3(z) - 4Li2(z) log z + 6C2 log z + - log^ z - 6(3 + 16Li2(z) - 8(z + 1)C2 



x1 e rr 32, 64(1 - z) 

81ogz-2(l-z)| + ^— ^|[-— logz+ ^ ^ 



log3(l - z) 



3 " ■ 3 

+ [ - 32Li2(z) - log2 z + 32C2 + 32(1 - z)] log2(l - z) + \- 40Li3(l - z) 
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+ 6Li4 ^1 - + 30Li4(l - z) + 26Li4(z) + 3Li2(z) log^ z - 10Us{l - z) log z 
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